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Abstract: 

It is proved that a source of electromagnetic radiation cannot emit a 
diffraction-free beam at the wave zone. A Bessel Jq 93-invariant beam does 
not hold even at the intermediate zone. These results negate claims published 
recently in the literature. 



An idea of creating a diffraction-free beam has been publishedfl]. The 
beam's amplitude is cylindrically symmetric ((^-invariant) where the r-dependence 
is proportional to the Bessel function of the first kind Jo(ar) and a is a con- 
stant having the dimension of L -1 . Reference [1] has arisen a great interest in 
utilization of Jo beams and has been cited more than 360 times[2]. An appli- 
cation of [1] shows the central peak of the assumed Jo beam [3] and another 
one refers to its peculiar ^-component wavelength [4]. Another publication 
related to [1] claims that a superluminal propagation of light in air has been 
detected [5]. Objections to [5] have been published[6]. The purpose of this 
work is to show that one cannot construct a diffract ion- free electromagnetic 
beam at the wave zone and that the Bessel function Jo(ar) is unsuitable 
for describing diffraction free (^-invariant wave at the intermediate zone too. 
This outcome proves that results of papers discussing this topic, in general, 
and those ascribing superluminal velocity to beams that take the form of 
Bessel function J , in particular, should be reevaluated. Units where the 
speed of light c = 1 are used. The metric g a p is diagonal and its entries are 
(1,-1,-1,-1). u r , rip and u 2 denote unit vectors in cylindrical coordinates and 
u x , u y and u 2 are unit vectors in Cartesian coordinates. 

A general analysis of diffraction-free solutions of Maxwell equations has 
been published [7]. Here the fields solving the problem are derived from a 
vector potential A that satisfies the wave equation together with the Lorentz- 
gauge requirement [8] . It turns out that this work is relevant to [1] and some 
of its results are analyzed here in detail. (Another work[9] is closely related 
to [1] and [7].) Let us start with the solution obtained in Example 1 (on 
p. 1557 of [7]). Using cylindrical coordinates and removing constant factors, 
the time dependent monochromatic electric field of this solution is obtained 
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from the vector potential E = —dA/dt 

E = uJx{ar)e< hz -^n v (1) 

where J\ is the Bessel function of the first kind of order 1. The magnetic 
field is B = curlA 

B = -bJ^ar) e^-^xir - iaJ Q {ar)e i{bz - ut) u z . (2) 

Ignoring constant factors, one finds that the magnetic field (^) is dual to 
the electric field of Example 2 of [7]. (The factor 2 in U r of example 2 is 
a misprint.) This outcome indicates that Examples 1 and 2 of [7] represent 
dual electromagnetic solutions where E — > B, B — > — E (see [8], p. 252). 

Having the solution, let us examine the problem of a cylindrically shaped 
wave guide whose walls are made of a perfect conductor (see [8], p. 335). 
The length of the cylinder is much greater than both its diameter 2R and the 
wavelength A = 1/uj (see fig. 1). The boundary conditions along the wave 
guide's walls are (see [8], p. 335) 

E { \ =0, B ± = 0. (3) 

Thus, the solution ([!]) and (0) satisfies the boundary conditions provided 

J^aR) = 0. (4) 

Dynamical properties of the solution ([!]) and ([|) are obtained from the 
energy-momentum tensor of the electromagnetic fields (see [10], p. 81 or [8], 
p. 605)) 

T f = ^(F» a F^g a p + ^FadsT) (5) 

where F^ u denotes the tensor of the electromagnetic fields. Expression (|^) 
is quadratic in the fields. Hence, one should use the real part of ([!]) and 
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(0) in an evaluation of quantities belonging to it. Let us first examine the 
momentum density of the fields. This is the Poynting vector 

S = — E x B. (6) 

The z-component of the momentum density and energy flux are obtained 
from the substitution of the appropriate real part of ([[]) and (fj) 

S z = —Ji(ar)cos 2 (bz — tot). (7) 

4:71 

Expression ([7]) is non-negative at all points, a property which is consistent 
with the beam's expected flux of energy that travels away from a localized 
source. 

The radial component of the momentum density is obtained analogously 

S r = — — Ji(ar)J (ar)sin[2(kz — oot)]. (8) 

Here one sees that, unlike the case of (0), the sign of (|8|) alternates periodically 
in the time and ^-coordinates. Moreover, for any fixed value of t and z, it 
changes sign along the r-axis, because zeroes of the Bessel functions J and 
J\ do not coincide [11]. It follows that although the radial motion does not 
vanish locally, its mean value is null. This property indicates that the radial 
motion takes the type of a standing wave. 

Now let us examine the interaction of the fields with the walls of the wave 
guide. Point P at x = R, y = z = is used as a representation of the general 
case and cartesian coordinates are used. The x-component of the momentum 
current at P is (see [10], p. 82 or [8], p. 605)) 

T xx = ^(E 2 y + El - El + Bl + Bl - Bl). (9) 

Examining the fields (]]]) and (0) and the boundary value (^), one finds that 
only the z-component of the magnetic field makes a nonvanishing contribu- 



tion. Thus, the momentum current at P is 



a 2 

T xx = — J 2 {aR)sin 2 {bz - ut). (10) 

This momentum current is absorbed by the walls, because the fields vanish 
in all space outside the inner part of the wave guide. 

Another effect of the magnetic field (0) on the wave guide's walls is the 
electric current induced in the ^-direction. Indeed, let us evaluate the line 
integral along the infinitesimal rectangular closed path of fig. 1. Using vector 
analysis, Maxwell equations and the boundary condition (Q), one finds 

<^B-dl = J cur IB- ds = J Air'yds. (11) 

Thus, a nonzero current j is induced on the walls, because only B z at the 
inner part of the wave guide makes a nonvanishing contribution to the line 
integral. This outcome proves that a time-dependent (and z-dependent) 
electric current flows along the (/9-direction of the wave guide's walls and that 
fields of this current are part of the solution (0) and @. This electric current 
sustains the B z related standing wave in the radial direction. The walls also 
counteract against local electromagnetic pressure. 

The dual solution of example 2 of [7] behaves analogously. Using the 
same global factor of ([!]) and @, one finds for this case 

B = uJ l {ar)e i ^ hz -^n tp (12) 

E = bJiiar^-^xir + iaJ {ar)e i{hz - wt) M z . (13) 
Hence, the boundary conditions (||) yield 

J (aR) = 0. (14) 

Since Jo{ar) and J\(ar) have no common root [11], a nonvanishing radial 
electric field exists at the wave guide's walls. It follows from Maxwell equation 



div~E = Anp that a time dependent and ^-dependent charge density is built 
on the inner part of the wave guide's walls. Thus, we have also in Example 2 
a current that flows on the walls and affects the fields inside the wave guide. 

Let us examine an analogous experimental setup. Here the source of the 
radiation at z — —L is the same as that of the first experiment but the 
wave guide is removed. This situation is different from the wave guide case. 
Indeed, the fields of a closed electromagnetic system depend on charges and 
currents at the retarded space-time points (see [10], pp. 158-160 or [8], p. 
225). Therefore, the wave guide's solutions clearly do not hold for this case 
because here the current along the wave guide walls is missing. 

Since in the second experiment the region at z = satisfies the wave zone 
requirements (see [10] p. 170 or [8], p. 392) 

L>A, L>2#, (15) 

one can use the wave zone solution. Let A denote the retarded vector po- 
tential at the wave zone. Thus, one finds the fields (see [10] p. 171) 

B = Axn, (16) 

E=(Axn)xn (17) 

where n is a unit vector in the radial direction. 

It turns out that the solution for the free space experiment is inherently 
different from the one which fits the wave guide's inner space. In particular, 
in the case of free space, fields at the wave zone are perpendicular to the 
radius vector from the source to the field point. On the other hand, the wave 
guide solution contains a ^-component (B z or E z ) which is an inherent part 
of the solution. As shown above, the B z (or E z ) field is associated with the 
electric current induced on the wave guide's walls. This conclusion obviously 

6 



holds for any pattern of source elements put at the same spatial region as 
the one used here, because the analysis does not refer to the source's details. 
Thus, the results disagree with the claim of [9]. 

One can use general arguments for proving that a diffraction-free electro- 
magnetic beam that has a nonvanishing ^-component for at least one of the 
fields, contains transverse standing wave. Indeed, the beam carries energy 
and therefore S of @ does not vanish. Hence, E is not parallel to B and, due 
to the z-component of the fields, S has a nonvanishing transverse component. 
Now, the diffraction-free property of the beam prevents energy from flowing 
transversally. Hence, the transverse component of S is a standing wave. 

It can also be proved that all solutions of [7] have a nonvanishing z- 
component of at least one of the fields. Indeed, the vector potential A takes 
the form (see p. 1556 therein) 

A = X)KM n + A»N n ), (18) 

n 

where a n and f3 n are numerical coefficients of the expansion. Here 

M n = curl[J n (ar)e t{bz+nip ~ wt) u z ] (19) 

and 

N n = jCurlM n (20) 

K 

where k is the wave number. Now N n contains a z-component (see p. 1557 
therein). Hence, if (3 n ^ then E = —dA/dt = iuA has a ^-component 
too. In other cases all (3 n = 0, which mean that for at least one n, a n ^ 0. 
Here the magnetic field B = curl A = a n curWl = fca n N, which means that 
B z ^ and the proof is completed. 

It follows that the family of solutions of [7] involves standing waves as- 
sociated with the components of the solutions. This diffraction-free family 
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of solutions may fit cylindrical wave guides but are unsuitable for the case of 
a free space. 

Example 4 of [7] (see p. 1558) is the last one which is analyzed here in 
detail. This example contains one component which is proportional to Jo(ar) 
and is yj-invariant. Although it has a ^-dependent z-component term which 
is associated with a standing wave, it looks simpler to show another problem 
of this solution. The vector potential of this example is given in Cartesian 
coordinates 

a 2 

A = -ia[aJ (ar) u x - i—J x {ar)cos^M^e< hz -^. (21) 
Using E = —dA/dt, one finds 

a 2 

E = auj[aJ (ar) u x - i—J x {ar)cos(p uje*^""*'. (22) 

Let us examine the ^-component of the Poynting vector which represents 
energy current flowing along the beam's direction, namely, the quantity which 
is analogous to ([?[) of Example 1 . Examining (|22| ) , one finds that only B y is 
needed for this purpose. Thus, (curl A) y of (|2~ID is 

3 3 

By = a[(ab-^)J (ar) + ^cos2<p J 2 (ar)]e« bz -^ . (23) 

Hence, the required ^-component of the Poynting vector is obtained as the 
product of the real parts of E x of (|22| ) and B y of (p3|) 

4 4 

S z = a 2 to[(a 2 b — -zT)Jo(ar) + —cos2ip J (ar)J2(ar)]cos 2 (bz — cut). (24) 

Let us examine the ^-component of the energy current near a point whose 
radial coordinate is R and Jq{clR) = 0. In this neighbourhood J 2 is domi- 
nant [12] and the contribution of the Jl(ar) term of ( |2"4| ) can be ignored. The 
rest of ( |24]) is proportional to jQ(ar)J 2 (ar)cos2ip. Now, let us examine the 
value of S z on a circle whose radius is R + s, where e is an appropriate small 



quantity. Due to the factor cos 2tp, one realizes that S z takes different signs 
on this circle. Hence, in the solution of Example 4 of [7], energy flows in 
opposite z-directions in certain regions of space. This property of Example 
4 is inconsistent with the notion of a beam, where electromagnetic energy 
flows away from a localized source. 

It is clear from the analysis carried out above that, in free space, one 
cannot build a diffraction free beam from the family of Bessel functions of 
[7], because these functions are unsuitable at the wave zone. 

Some conclusions can be drawn for the intermediate zone too. The diffrac- 
tion free (^-invariant Jo(ar) function proposed in [1] does not belong to the 
solutions of [7]. Indeed, in [7], there are only two truly (/2-invariant solutions. 
They are the dual solutions of Examples 1 and 2 which are discussed above. 
As proved in this work, the ^-component of the energy current is propor- 
tional to Jf(ar). Hence, the flow of energy vanishes along the z-axis. It is 
also proved above that Example 4 of [7], where there is one J term which 
is (/9-invariant, does not describe a beam of electromagnetic radiation and 
its ^-component is not (/9-invariant. It follows that experiments using a tp- 
invariant setup and showing a strong peak at the center (like [1,3,4]) should 
not be interpreted by means of diffraction free solutions. 
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Figure captions: 



Fig. 1: 

Electromagnetic radiation is emitted from a source into a cylindrical wave 
guide whose radius is R. The source is at z = —L and L S> 2R. O denotes the 
origin of coordinates and the rectangle at point P denotes a closed integration 
path (see text). 
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